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ABSTRACT 

The idea that cosmic relativistic jets are magnetically driven Poynting-dominated flows 
has many attractive features but also some problems. One of them is the low efficiency 
of shock dissipation in highly magnetized plasma. Indeed, the observations of gamma 
ray bursts (GRBs) and their afterglow emission indicate very high radiative efficiency 
of relativistic jets associated with these phenomena. We have revisited the issue of 
shock dissipation and emission and its implications for the internal shock model of the 
prompt GRB emission and studied it in the context of impulsive Poynting-dominated 
flows. Our results show that unless the magnetization of GRB jets is extremely high, 
a > 100 in the prompt emission zone, the magnetic model may still be compatible 
with the observations. First, for a ~ 1 the dissipation efficiency of fast magnetosonic 
shock is still quite high, ~ 30%. Second, the main effect of reduced dissipation effi- 
ciency is merely an increase in the size of the dissipation zone and even for highly 
magnetised GRB jets this size may remain below the external shock radius, provided 
the central engine can emit magnetic shells on the time scale well below the typical 
observed variability scale of one second. Our analytical and numerical results suggest 
that magnetic shells begin strongly interact with each other well before they reach the 
coasting radius. As the result, the impulsive jet in the dissipation zone is best described 
not as a collection of shells but as a continuous highly magnetised flow with a high 
amplitude magnetosonic wave component. How exactly the dissipated wave energy is 
distributed between the radiation and the bulk kinetic energy of radial jets depends 
on the relative rates of radiative and adiabatic cooling. In the fast radiative cooling 
regime, the corresponding radiative efficiency can be as high as the wave contribution 
to their energy budget, independently of the magnetization. Moreover, after leaving 
the zone of prompt emission the jet may still remain Poynting-dominated, leading to 
weaker emission from the reverse shock compared to non-magnetic models. Energeti- 
cally sub-dominant weakly magnetized "clouds" in otherwise strongly magnetised jets 
may significantly increase the overall efficiency of the shock dissipation. 

Key words: MHD — relativity — gamma-rays: bursts — ISM: jets and outflows — 
galaxies: jets 



1 INTRODUCTION 



of the reasons why this mechanism has attracted so much 
attention in recent years. 



Various astronomical observation indicate, directly or indi- 
rectly, the existence of highly relativistic outflows in a va- 
riety of cosmic phenomena, such as active galaxies, pulsar 
wind nebulae, X-ray binary stars, and Gamma Ray Bursts 
(GRBs). Although the origin of these flows is still a subject 
of debate, especially in the case of GRBs, and requires fur- 
ther investigation, so far we have identified only one mech- 
anism of jet production which may operate in all these very 
diverse environments - the magnetic mechanism. This is one 



It has been shown that magnetic fields can not only 
tap the rotational energy of a massive rotator, placed in the 
"hart" of "cosmic jet engines" in this model, but also ac- 
celerate and collimate outflows. In fact, it has been shown 
that in the case of relativistic steady-state jets, their mag- 
netic acceleration and collimation go hand-in-hand and for 
this reason this version of magnetic mechanism is called the 
collimation acceleration mechanism. In order to ensure effi- 
cient magnetic acceleration, the jet opening angle should be 
small compared to the jet Mach angle associated with fast 
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magnetosonic waves (jKomissarov et al"1l2009ah . In the small 
angle approximation, the Mach angle 

$M — ~r~7 i 

M 

where M is the fast magnetosonic Mach number. In the case 
of predominantly asymuthal magnetic field this number can 
be estimated using Eqf3]and the acceleration condition reads 

° > irA)\ 

where a — B 2 / pc? . In order to understand where this lim- 
itation comes from, consider a freely expanding conical jet. 
If its magnetic field is predominantly azimuthal than the 
magnetic freezing yields B a r _1 and the magnetic energy 
£ m oc r 2 B 2 oc r° is not utilised to accelerate the flow. The 
flow geometry has to deviate from the conical one in order for 
the magnetic acceleration to operate and this requires effi- 
cient causal communication across the flow fe.g. iKomissarovl 
l201ll ). In fact, the above causality condition is a bit too strict 
and the acceleration may proceed, though at a much lower 
logarithmic rate, even after it is no longer satisfied. How- 
ever, t he above constrain t on a remains valid up to a factor 
of few (|Lvubars kv 2009). For AGN jets with their inferred 
7j ~ 10 and observed 0j ~ 0.1 this yields a > 1, whereas for 
the commonly accepted parameters of GRB jets, jj > 100 
and 9j ~ 0.1, the last equations implies a 3> 1. 

The observed variability of the emission produced in 
these jets has been traditionally associated with strong 
shock waves, driven into the jets by the ir unsteady cen- 
tral engines (see the review bv lPiranlBobl ). Their observed 



bright knots and spots have also been attributed to such in- 
ternal shocks. Indeed, shocks are generally known as places 
of effective dissipation of kinetic energy and acceleration of 
relativistic electrons, responsible for the non-thermal emis- 
sion observed in many astrophysical objects. However, in 
the case of relativistic magnetized flows this interpretation 
encounters sig nificant problems, pa rticularly in the case of 
GRB jets (e.g. iNaravan et alj|201lh . 

Indeed, the dissipation efficiency of relativistic shocks 
in highly magnetized plasma is rather low. This was 
widely accepted already after the pioneering work by 
iKennel fc CoronitH (1984). To be more accurate, this state- 
ment is concerned with fast magnetosonic shocks only. Slow 
magnetosonic shocks can still have very high dissipation ef- 
ficiency, because at such shocks the magnetic energy can be 
dissipated as well l|Lvubarsky||2005l ). However, fast shocks 
are much more readily produced, usually via collisions, 
whereas formation of slow shock requires some rather special 
conditions. 

In contrast, the observations indicate that the jet ra- 
diative efficiency, rj r , defined as the fraction of jet en- 
ergy eventually converted into radiation (usually non- 
thermal), can be quite high. For example, the observa- 
tions of the GRB afterglows imply that the radiative ef- 
ficiency of their jets is often i n excess of 10% and some- 



times may even reach 90% JPanaitescu Sz Kuniai 20021 ; 



ay 

lYost et all 120031 ; iGranot et all l200d ; IZhang et all I2007al 1. 

According to the more recent study of Swift GRBs the 
situation is even more dramatic, with the me an radiative 
efficiency around 90% (|Willingale et all 120071 ). This diffi- 
culty of the shock model has forced many theorists to start 
looking for alternative models involving direct dissipation 



of magnet ic energy associated with t h e magnetic reconnec- 
tion fe.g. iDrenkhahn fc Spruij |2002|; IZhang fc Yan I I20T1I ; 
iMcKinnev fc Uzdenskvl l20ld ; iLvubarskvl |2010l )." However, 
these models still remain at a rather rudimentary level of 
development due to difficulties of their own. 

In most of the previous theoretical studies of 
magnetized relativistic jets it was assumed that 
they were more or less homogeneous, just for sake 
of simplicity. The strong observed variability and 
the complex observed structure of some relativis- 
tic jets, where high resolution images are avail- 
able, suggest that this may be not a very re- 
alistic assumption. Following the early work by 
IContopoulosI (|l995T ). a number of recent papers ex- 
plored the implications of hig hly intermittent jet 
production on its dynamics ([Gr anot et al. 2010; 
Lvutikovl l2010l: iLvutikov &; Lister! |2010| : ILvutikovl 



2011]; 



Granot 



2011c 



They have concluded that 
longitudinal expansion of highly magnetised plasma 
shells may result in efficient conversion of the Poynt- 
ing flux into the bulk kinetic energy of the shells and 
strong reduction of their magnetisation. However, in 
those papers only the dynamics of a single shell was 
studied in details, whereas the case of an impulsive 
jet composed of many such shells was subjected to 
a much more speculative analysis. The main goal of 
our study was to reduce this imbalance. 

As a first step in studying of the multiple shell 
case one may assume that the gaps between them 
are empty. However, if the jet engine does indeed op- 
erate in an impulsive fashion then external plasma, 
presumably of much lower magnetization, is likely 
to fill the gaps during quiescent periods. Then each 
time the jet is reborn it has to push this plasma 
aside. Provided the jet is sufficiently powerful, this 
can be done quite efficiently b y the bow-shock de- 
velop ing at the jet head (e.g. iKomissarov &; Fallel 
1998). Farther out, where the distance between 
shells becomes comparable to the jet radius, the 
relativistic effects make impossibl e for the externa l 
plasma to enter the gaps (|Lvutikov & Lister! IbplOh . 
At the same time it becomes impossible for the en- 
trained plasma to leave the gaps as it is forced to 
remain within the cone of the half-opening angle 
1/Yj. Thus, it seems quite plausible that some of the 
external plasma will remain in the gaps and become 
part of the jet, though it is still rather difficult to 
quantify the effectiveness of this mass-loading pro- 
cess at present. 

On one hand, loading Poynting-dominated jets 
with weakly magnetised clouds complicates the 
problem. On the other hand, this may play a very 
important role in their physics. As far as the ra- 
diative efficiency is concerned, these clouds could 
be the locations there most of the shock dissipa- 
tion and emission takes place. Very much in the 
same way as in the model of the afterglow emis- 
sion, where a magnetic piston drives the so-called 
external shock wave th rough weakly magnetised in- 
terstellar medium (e . g. ILvutikov &z Blandfordl 2003; 
IZhang 8c Kobavashil 120051 : ILvutikovl l2006h . Obvi- 
ously, these clouds may have to be "excited" many 
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times before a significant fraction of magnetic en- 
ergy is radiated. 

A si milar repetitive "pumping " action has been investi- 
gated bvlKobavashi fc Saril ([2001) in the case of unmagne- 
tized highly variable jets. As individual shells (portions of 
the jet moving with very different Lorentz factors) collide 
and heat-up, only a fraction of the dissipated energy is radi- 
ated. The rest of it remains initially in the form of heat, but 
later, when the shells begin to expand, this heat is converted 
back into the kinetic energy of relative motion. When an- 
other collision occurs, a fraction of this energy is dissipated 
and radiated again and so o n. The process continues until 
the shocks become very weak. Kobavas hi fc Saril (|200ll ) have 
demonstrated that this way the radiative efficiency can be 
increased up to 60%, even if during each individual collision 
only 10% of the dissipated energy is radiated. However in 
order to achieve this, they required very strong variations of 
the jet Lorentz factor, with uniform distribution of log 10 7 
between 1 and 4. Since in our case the energy behind this dy- 
namics is of magnetic nature, the name "magnetic pump" 
springs to mind. In fact, each time the shock-heated gap 
plasma expands and its components, which are unable to 
cool radiatively, cool adiabatically their thermal energy is 
returned back to the "pump" and recycled. 

The whole problem of impulsive jet dynamics from its 
production to its interaction with the interstellar matter is 
still prohibitively complex. In order to make progress, we 
will consider much simpler problems hoping to elucidate 
some of its important aspects. We start with the issue of the 
dissipation efficiency of fast magnetosonic shocks in highly 
magnetized plasma as we feel need to clarify few important 
points. This is done in Section [2] and in Appendix [Al Then 
we consider strictly periodic one-dimensional flows in slab 
geometry in the framework of one fluid MHD with simple 
polytropic equation of state. In Section [3] we study oscil- 
lations developing in a periodic train of initially stationary 
magnetic shells, as they expand, collide and radiatively cool. 
Then we consider mov i ng tra ins with initially empty gaps, 
following iGranot et al.l f2010), first in the adiabatic regime 
(Section 0J) and then in the regime of fast radiative cool- 
ing (Section [5]). Finally, we study the dynamics of a moving 
train with gaps filled with weakly magnetised plasma from 
the start in the fast radiative cooling regime (Section^. The 
results and their astrophysical implications are discussed in 
Section [7] Our conclusions are listed in Section [5] 

Throughout most of the paper we use the Heaviside 
units, where c = 1 and B/\/4^r —¥ B, but in the Discussion 
we reintroduce the speed of light. 



2 DISSIPATION EFFICIENCY OF 

PERPENDICULAR FAST MAGNETOSONIC 
SHOCKS. 

Here we consider only the perpendicular relativistic shocks, 
where the flow velocity is perpendicular to the shock front. 
In addition we assume that the magnetic field is parallel 
to the shock front. Thus, we restrict our attention to purely 
one-dimensional flow directed perpendicular to the magnetic 
field. This constraint prohibits slows magnetosonic waves 
of any kind and the only non-trivial shocks solutions are 
the fast magnetosonic ones. Such shocks were first studied 



bv lfCennel fc Coronitil (|l984l ) in application to the termina- 
tion shocks of pulsar winds. In order to simplify the shock 
equations, they only considered the case of cold upstream 
flow, ultrarelativistic shock speed, and high shock strength, 
in the sense that the downstream Lorentz factor is much 
lower compared to the upstream one. All these additional 
constrains are justified in the case of pulsar wind nebulae. 
The general case of mag netos onic shocks was anal y sed b y 
iMaiorana fc Anild dl987h a n d lApple fc Camenzindl ljl988l l. 
La ter. [ Zhang fc Kobavashil ( 12005:) expanded the analysis 
of iKennel fc Coronitil (|l984T ) by allowing variable shock 
strength, which they described by the Lorentz factor of 
the relative motion between the upstream and downstream 
states, 712. Although this parameter can indeed be used to 
describe the shock strength, the more traditional parameter, 
unanimously accepted in the non-relativistic hydrodynamics 
and MHD, is the shock Mach number. The proper relativis- 
tic definition of Mach number with respect to the wave mode 
of speed c m in the fluid frame is 



M = 



7 ' 1 



where 7 m is the Lorentz factor corresponding to the wave 
speed c m and 7 is the Loren tz factor corresponding to the 
flow speed v (iKonigl 1 1 1980T ) . In the limit of cold plasma, 
where the thermodynamic pressure, p, and hence the sound 
speed are set to zero, the fast magnetosonic speed is the 
same in all directions 

B 2 



2 

c f = 



B 2 + p 1 + a' 



(2) 



where B and p are the magnetic field strength and the gas 
rest mass density as measured in the fluid frame, and 

* = - (3) 
P 

is one of the parameters describing the plasma magnetiza- 
tion. For 7 1 and a 3> 1 the fast magnetosonic Mach 
number is 



M 



7 



(4) 



In Appendix |A"| we redo the analysi s of perpendicular fast 
shocks of lZhang fc Kobavashil (120051 ) using the Mach num- 
ber with respect to the fast magnetosonic mode as the shock 
strength parameter. While in general numerical techniques 
have to be used to solve the shock equations, in the limit of 
high shock Mach number, Mi 3> 1, and high upstream mag- 
netisation, <ti ^> 1, they allow simple approximate solution 
where 



72 
pi 

P2 



1/2 
Mi/91, 

\piMl 



(5) 
(6) 

(7) 

(8) 



where p m is the magnetic pressure. Here index "1" refers 
to the upstream and index "2" to the downstream state. 
One can see that for Mi 3> 1 there are strong jumps in the 
rest mass density and magnetic pressur e as measured in the 
fluid f rame. This is what is meant by IZhang fc Kobavashil 
(2005), when they state that high magnetisation does not 
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Figure 1. Left panel: Dissipation efficiency of perpendicular shock as a function of the magnetisation parameter <r of the upstream state 
for three different values of the shock fast magnetosonic Mach number. The solid lines show the solutions for the polytropic equation of 
state with T = 4/3 (k = 4). The dashed lines show the solutions for the electron-positron Synge gas. Right panel: Dissipation efficiency 
of perpendicular shock as a function of the shock fast magnetosonic Mach number Mi for three different values of the magnetisation 
parameter a of the upstream state. The equation of state describes polytropic gas with T = 4/3 (k = 4). 



prevent development of strong shocks. On the other hand, if 
we consider parameters measured in the shock frame, which 
will be indicated by prime, then 



P2 r - 

B' 2 



Pi , 



(9) 
(10) 



where p' is the rest mass, not the inertial mass, density. Eq[9] 
shows that there is no much decrease in the shock frame 
volume occupied by plasma as it crosses the shock - the 
large decrease in the proper specific volume is almost totally 
compensated by the reduced Lorentz contraction. It is this 
what is meant when shocks in highly magnetised medium 
are often described as weak or weakly compressive. 

One can define the shock dissipation efficiency in many 
different ways, some more meaningful than others. We are 
interested in the fraction of the total energy flux which 
can be converted into radiation without invoking any ad- 
ditional dissipations mechanisms, like the magnetic recon- 
nection, downstream of the shock. This suggests to define 
the efficiency as 

F t ,2 



1b 



Jtot,2 



(11) 



where Ft = e t j 2 v is the thermal energy flux density and 
-ftot = v is the total energy flux density, e t is 

the thermal energy density and w = p + e t is the relativis- 
tic enthalpy, both defined in the fluid frame. Notice that 
this definition makes ?7s independent on the flow velocity, 
and hence on the shock speed relative to the observer. For 
Mi , <Ti S> 1 we find 



1 



2(1 + <Tl) 

(see Appendix 1X1) 



(12) 



Equation (|12[) shows that for o\ 2> 1 only a rather small 
fraction of the flow energy can be dissipated and then ra- 
diated away. In order to see if the same conclusion applies 
to shocks in only moderately magnetized plasma we solved 
the shock equations numerically. The results are presented 
in Figures {TJ and ([2]). The left panel of Fig.Q shows the 
dissipation efficiency as a function of <7i for three differ- 
ent values of the shock fast magnetosonic Mach number, 
Mi = 2, 10, and 100. One can see that for o\ 1 the dis- 
sipation efficiency does indeed decline as . However, for 
small magnetization the efficiency actually increases with o\ . 
The location of the maximum depends on the shock Mach 
number and for weak shocks is near o\ = 1. However, its 
magnitude is rather low in this case. As one can see in the 
right panel of Fig[T] the efficiency monotonically increases 
with the shock Mach number. 

The collimation acceleration of magnetic jets may result 
in the asymptotic magnetization a ~ 1. For this reason we 
presented in the right panel of Fig[T]the dissipation efficiency 
as a function of the shock Mach number for o\ around this 
value. For very high Mach number the efficiency is rj s ~ 0.8, 
0.3, and 0.05 for a — 0.1, 1, and 10 respectively. This shows 
that for <7i ~ 1 the dissipation efficiency can be already 
reasonably high. Moreover, Eq |12l gives a rather accurate 
estimate of the efficiency for Mi > 10 and <Ti > 1. 

The shock solution depends on the plasma equation of 
state (EOS) but not strongly, at least for the explored range 
of parameters. In the right panel of Fig[T] the solid lines 
show the solution for the polytropic equation of state with 
r = 4/3 (k = 4), whereas the dashed lines show the so- 
lution for the electron-positron Synge gas, which assumes 
the same tem perature rela tivistic Maxwell distribution for 
every species (|Svnge|[l957l ). One can see that difference is 
relatively small, particularly for strong shocks. We also ex- 
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Figure 2. Fractions of the upstream kinetic energy converted at the shock into the thermal energy, ft, magnetic energy, / m , and 
remaining in the kinetic form, f^. The left and the middle panel show the fractions as functions of the upstream magnetisation for 
Mi = 10 and Mi = 100 respectively. The right panel shows them as functions of the shock fast magnetosonic Mach number for the 
upstream magnetisation <ji = 1. 



perimented with the electron-proton Synge gas and found 
that the solution was even closer to that with the polytropic 
EOS. 

The result (|12|) has a straightforward interpretation. 
First, only the kinetic energy dissipates at the shock. Second, 
the kinetic energy flux makes only 1/(1 + <7i) of the total 
upstream energy flux. Finally, for high Mi approximately 
one half of the kinetic energy dissipates into heat, and ap- 
proximately one half is converted into Poynting flux. This 
is illustrated in Fig[2]which shows fractions of the upstream 
kinetic energy converted at the shock into the thermal en- 
ergy, ft, magnetic energy, / m , and remaining in the kinetic 
form, / k . 



3 CHAMBER OSCILLATIONS 

Collisions between shells will produce reflected shock waves, 
similar to those created via shock reflection off a conducting 
wall. So the problem of shell interaction appears analogous 
to that of a shock bouncing off the ends of a closed tube 
that contains plasma with very diverse magnetization. Each 
time it crosses the low magnetisation domain a fraction of its 
energy is dissipated and radiated away, so the shock weakens 
bit by bit. 

Consider a one-dimensional flow confined within a 
chamber of length I. Suppose that initially the chamber is 
divided into two sections, of lengths l p and l g = I — l p . The 
first section, which we will call the "pulse" , is filled with uni- 
form highly magnetized cold plasma and the second section, 
which will be referred to as the "gap", is uniformly filled 
with plasma of lower magnetization and weaker magnetic 
field. When comparing this configuration with inhomoge- 
neous relativistic jet, one is tempted to identify the length 
I with the separation between two neighbouring shells as 
measured in the jet frame. 

In the purely electromagnetic version of this problem 
the gaps are empty and the pulses have a uniform distribu- 
tion of magnetic field with vanishing electric field. The solu- 
tion to this problem involves two identical electromagnetic 
pulses bouncing between the perfectly conducting walls of 



the chamber without decay. When the plasma magnetization 
is high, a > 1, we expect the MHD solution to be close to 
the electromagnetic one. However, the shock dissipation will 
gradually damp these oscillations. If the radiative cooling of 
the chamber plasma is indeed very efficient, it eventually 
relaxes to an equilibrium with uniform magnetic field and 
negligibly small temperature. This allows us to compute the 
total loss of energy from the system, and hence its radiative 
efficiency. 

3.1 Asymptotic state 

Denote as p p , B p , and M p the rest mass density, the mag- 
netic field, and the total mass of the pulse respectively, and 
as p g , B g , and M g the corresponding parameters of the gap. 
It is convenient to describe the problem by the ratios of 
lengths, rest masses, and magnetic energies of the gap and 
the pulse: 



X _ h X — M 3 



g = B 2 g lg 

Bplp 



(13) 



The relativistic magnetization parameter of cold 
plasma, a = B 2 / p, gives the ratio of the magnetic and rest 
mass energies in the fluid frame. Given this, it makes sense 
to define the mean magnetization of plasma in the initial 
state as 



Mo 



1 B p lp+Bglg <J p + (TgS„ 



P P 1 P + pgig 



1 + S„ 



(14) 



The condition of strong mean magnetization constrains the 
mass fraction of the system. In particular, if the gap mag- 
netization is really low, a g <C 1, this requires a p ^> 8 m - 

From the rest mass conservation and the magnetic field 
freezing we have 



Pglg — Pgig i p P l P — p P l P 



Bglg — Big , 



B p l p — Bl p , 



(15) 
(16) 



where tilde denotes parameters of the equilibrium state, 
which is reached asymptotically for t — > oo. These combine 
to yeild 
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Figure 3. Radiative efficiency in the chamber problem for <5; = 
0.5 (solid line), <5; = 1 (short dash line), 5; = 2 (dot line), <5; = 3 
(dash-dot line), <5; = 9 (long dash line) as a function of the energy 
ratio parameter S e . 



B = B; 



(17) 



Using this result one can derive the magnetic energy of the 
equilibrium state and hence the radiative efficiency, 



{l + 5] /2 5V 2 Y 

(l+*)(l + *e) ' 



1 - 



(18) 



which is defined here as the fraction of the initial magnetic 
energy converted into radiation. This function is shown in 
Figure [3] One can see that the radiative efficiency increases 
with 8i and decreases with <5 e , which has a very simple ex- 
planation. Smaller <5 e means stronger expansion of the pulse 
and hence smaller magnetic energy remaining in the system 
after its relaxation. For 8 e — the efficiency is simply the 
fraction of the volume available for the pulse to fill, 

»7r,max = y . (19) 

The gap plasma however resists the pulse expansion and 
its resistance increases with the gap pressure and hence the 
magnetic energy stored in the gap. 

This simple analysis hints that the radiative efficiency 
of impulsive magnetically dominated flows can be very high. 
Moreover, the outcome does not even depend on the magne- 
tization. Other important properties of the process however 
may do. For example the rate of dissipation, and hence the 
luminosity. In order to investigate this issue a bit further we 
have carried out numerical simulations. 



3.2 Numerical simulations 

We solve numerically the one dimensional equations of 
single-fluid relativistic magnetohydrodynamics in plane (or 
slab) geometry. In order to account for the radiative cool- 
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Table 1. Parameters of numerical models for the chamber prob- 
lem. (<r)f is the final overall magnetization of plasma, r; r , max is 
the maximum radiation efficiency given by Equation 1181 and tso 
is the time by which the plasma emits 50% of the energy corre- 
sponding to this efficiency and given in the units of the chamber 
light crossing time. Other parameters are explained in Section [3] 



ing we introduce a source term in the energy-momentum 
equation: 



qu 



(20) 



where T" M is the stress-energy-momentum tensor of fluid, 
u M is its 4-velocity, and q is the cooling rate as measured 
in the fluid frame. The cooling rate used in our simulations 
is not based on any particular physical mechanism. At this 
stage, it makes sense only to require for the cooling time 
to be short compared to the dynamical time scale, as this 
seems to be required by the observations of GRBs. To be 
specific, we put 



fc(T) 



Ar c , 



(21) 



where T = p/pc 2 is the "temperature" 



0, 

(T-ToJ/To, 
1, 



T < To 

T < T < 2T 

T > 2T 



and At coo i is the characteristic cooling time. In the simula- 
tions we used To = 10 -3 and At CO oi = 0.04(£/c), a small 
fraction of the chamber light crossing time. This cooling 
function implies that all particles cool rapidly and this may 
be rather unrealistic, as only relativistic electrons radiate 
efficiently. However, we do not expect this to change the 
results by more th an a factor of few becau se the 'thermal 
pump' mechanism l|Kobavashi fc SarH [2001) is likely to be 
efficient in re-processing the retained thermal energy. The 
equation of state describes an ideal gas with the adiabatic 
index T = 4/3. 

The nume r ical s cheme is based on the one described 
in iKomissarovl (|l999T ) , with few improvements added over 
the years. The computational grid is uniform with 600 cells. 
Because of the presence of strong shocks in the solution the 
scheme is only of the first order accuracy. Our convergence 
study shows that the typical computational error for most 
of the parameters discussed below is about few percents. 
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Figure 4. Magnetically driven oscillations in the model X of the chamber problem. From left to right, the plots show the solution at 
t = 0.2, 0.42, 0.9 and 1.5 respectively. The first row shows the magnetization parameter <r, the second row the magnetic pressure p m , 
the third row the flow velocity, u x = -fv x , and the bottom row shows the gas temperature T = p/p. 



The units are selected in such a way that the dimen- 
sionless chamber length is I = 1 and the speed of light is 
c — 1. The pulse is initially located in the middle of the 
chamber. The transition from the pulse to the gap state is 
smoothed out using the tanh(:r) function, with the width of 
the transition layer 8x — 0.03. 

Table Q] shows the parameters of all models we present 
here. Following the reasons behind the magnetic pump 
mechanism, the magnetization parameter of the pulse, a p , is 
always high, varying between 5 and 30. The gap magnetiza- 
tion is normally much lower, a g ~ 0.1, but we also included 
the cases where a g is above unity. In all cases, the total en- 



ergy of the initial state is dominated by the magnetic field, 
corresponding to (<r)o > 1. 

We first discuss the model X, where the gap is basically 
e mpty, as this is the c losest case to the scenario envisioned 
in lGranot et alj (|2010T ). Figure [3] illustrates the flow dynam- 
ics in this model. The first column shows the solution at 
t — 0.2, when the pulse is beginning to expand for the first 
time. In both halves of the chamber, the solution exhibits 
a rarefaction wave, moving into the pulse and ejecting its 
plasma into the gap. The ejected plasma expands into the 
gap space and develops very high velocity. In fa ct, it can 
reach 7 = 2ctq + 1 when t he gap is a pure vacuum (|Lvutikovl 
|2010| ; iGranot et aLll2010l ) , but in our simulation it is limited 
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by the numerical resolution and the non-vanishing density 
of the gap plasma. 

The second column shows the solution at t = 0.42, soon 
after the collision of the flow with the chamber walls, which 
drives strong shock waves back into the pulse. One can see 
that behind the shocks both the magnetic pressure and the 
plasma magnetization are almost as high as initially in the 
pulse. The gas passed through the shock is heated to a very 
high temperature and its cooling rate reaches maximum. On 
the contrary, in the initial location of the pulse a secondary 



gap is forming behind the two rarefaction waves which have 
now been reflected off the centre. The gap and the pulse have 
switched their places almost recreating the initial conditions. 

The third column shows the solution at t = 0.9, soon 
after the reflected shocks have collided in the center of the 
chamber for the first time. One can see that a region of high 
magnetic pressure and a is formed in the original location 
of the pulse. This "born again" pulse drives the next round 
of expansion and shock heating. 
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Figure 6. The rate of radiation losses normalized to Eo/tg, where Eq is the total initial magnetic energy and to = l/c is light crossing 
time of the chamber. 



The forth column of Figure ((4} shows the solution at 
t = 1.2, which is separated from the solution presented in 
the first column by exactly one light crossing time of the 
chamber. Qualitatively, the solutions look very similar, but 
at t — 1.2 the gap is no longer empty and the pulse expansion 
drives a shock through its plasma. This is why the Lorentz 
factor is so much lower. 

Figure [5] shows for comparison the solution for the 
model D, where the gap is filled with a significant amount 
of plasma from the start. One can see that there is a strong 
similarity between models D and X. The first column shows 



the solution at t = 0.2, when the pulse is beginning to ex- 
pand and compress the gap plasma for the first time. In 
both halves of the chamber the solution exhibits a rarefac- 
tion wave, moving into the pulse, and a shock wave, moving 
into the gap. At the shock, the plasma is heated and then 
it rapidly cools. An insignificant heating is also seen in the 
middle of the rarefaction wave, where the velocity gradient 
is the highest. This heating is entirely due to the numeri- 
cal viscosity, and represents computational errors. The sec- 
ond column shows the solution at t = 0.42, soon after the 
shocks reflection off the chamber walls. By this point the 
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Figure 7. Total energy radiated by the time t normalized to the 
initial magnetic energy in the chamber problem. The dashed line 
shows the model X, which initially has an "empty" gap. 

gap plasma has been shocked twice, first by the incident and 
then by the reflected shock. It's temperature and the cooling 
rate are reaching maximum. By t = 0.9, the reflected shocks 
have collided in the center of the chamber for the first time 
and the "born again" pulse drives the next round of expan- 
sion and shock heating. When the shocks reflect of the walls 
again (see the fourth column of Figure[5]l the solution is very 
similar to that at the time of the first reflection. Thus, the 
system is undergoing strong oscillations and the gap plasma 
experiences repeated shock-heating. 

Figure [6] shows the integral rate of radiative cooling for 
all our models, except G. These curves exhibit strong peaks 
at the time of the shock collisions at the walls and in the 
middle of the chamber, as this is the time when the shocks 
cross plasma with lower magnetization. In model X, the am- 
plitude of the peaks decreases with time monotonically. This 
reflects the fact that in this models the gaps developing in 
the centre of the chamber have very similar parameters to 
those developing at the walls. The overall decrease of the 
cooling rate is caused by the gradual decrease of the shock 
strength and increase of the gap magnetisation. In all other 
models, peaks corresponding to the shock collision with the 
wall are significantly stronger than those corresponding to 
the shock collisions in the middle of the chamber. This is 
because the gaps are filled with weakly magnetised plasma 
from the start. 

These peaks are strongest for models C,D,E and G, 
which have lower gap magnetization, a g = 0.1, and much 
weaker in models A,B,F, and H which have higher gap mag- 
netization, a g = 1,3. The integral radiative cooling rate of 
models with higher gap magnetization is also much lower 
and it declines slower. This is fully consistent with lower en- 
ergy dissipation rate in the case of higher magnetization, as 
explained in Section [21 and with the fact that the overall 
radiative efficiency does not depend on the magnetization 
(see Section [3]). 

Figure [7] shows the total energy SE(t)/Eo radiated by 
the time t. It confirms that even models with high gap mag- 
netization eventually approach the efficiency i) r . m „ given by 



Eq |18l this only takes a bit longer. The last column of Ta- 
ble [1] gives the time required for the system to radiate 50% 
of ^r.maxi^o, the total amount of energy which will be even- 
tually lost. It spans from one to ten chamber light crossing 
times, increasing as expected with the gap magnetization. 
It is interesting that the model X is not much different from 
models, in spite of the lack of weakly magnetized plasma in 
the initial state. The explanation for this seems to reside in 
the ability of strong rarefaction waves to supply such plasma 
and provide it with much higher fracti on of the kinetic en - 
ergy compared to the rest of the flow l|Granot et alj|2010h . 
This agrees with the fact that faster dissipation is found 
in the models C and G, which have three times wider ini- 
tial gaps, thus allowing a larger fraction of the ejected pulse 
plasma to develop low a and high 7. The model E also shows 
faster dissipation. This is due to the fact that its initial gap, 
which is already weakly magnetized (a = 0.1), is ten times 
more massive compared to other models. 



4 ADIABATIC FLOW 

Although suggesting many interesting hints, the chamber 
problem may differ from the case of impulsive flow in many 
respects. For example, in this problem all of the released 
magnetic energy is converted into radiation. This can not 
be so in the case of a flow, where not only energy but also 
momentum is associated with the electromagnetic field of 
pulses. Some of this momentum can be passed on to plasma, 
resulting in its bul k acceleration. Consi dering the dynamics 
of an isolated pulse. iGranot et al.l (|2010T ) have found that the 
pulse develops a very fast "head" and that most of the pulse 
magnetic energy is converted into the bulk motion kinetic 
energy of this head. In principle, interaction b etween pulses 
may c hange this outcome. Analysing this issue. lGranot et aU 
(2010) pointed out that the head does not spread much un- 
til its acceleration is over and argued that because of this 
the collisional interaction between individual pulses is de- 
layed until the end of their acceleration phase. However, each 
pulse also develops a long slow tail and even if the gaps be- 
tween pulses are completely evacuated initially they quickly 
become filled with the plasma of these tails. The shock in- 
teraction between the heads and the tails may modify the 
flow dynamics. 

The easiest way to address this issue is to consider a pe- 
riodic train of identical traveling pulses separated by empty 
gaps. As the pulses interact and the shocks dissipate the 
kinetic energy of relative motion, the flow gradually ap- 
proaches a state with uniform total pressure. Provided the 
radiative cooling is negligible, the parameters of this state 
can be found analytically from the equations of mass, mag- 
netic flux, energy, and momentum conservation. 

4.1 Asymptotic state 

These conservation laws are 

[pjl] = , (22) 

[Bl] = , (23) 

[( W7 2 -p + J B 2 (l-l/27 2 ))*] =0, (24) 

[(™ 7 2 + B 2 )vl] = , (25) 
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Figure 8. Various stages in the evolution of the adiabatic flow. The first column panels illu strate the solution at t = 0.5, when the 
pulse tail has just reached its head. This solution is similar to the solution for an isolated pulse iGranot et al.|[2010T) . The second column 
shows the solution at t = 1.0. At this point the double shock structure, which has been created by the collision of the head and the tail, 
is already beginning to have a strong impact on the flow. The third column shows the solution at t = 10. Now the flow is very different 
from the isolated pulse solution. The features of strong reverse rarefaction wave have been totally erased. The last column shows the 
solution at t = 1000. By this time the flow has developed the "saw-tooth" profile characteristic of the classic non-linear wave steepening 
problem. In each column the top plot shows the Lorentz factor, the middle plot shows the total pressure, which is dominated by the 
magnetic pressure, and the bottom plot shows the magnetization parameter a. 



where w = p + up is the relativistic enthalpy, k = T/(T — 1), 
where F is the adiabatic index, B and v are the magnetic 
field and the plasma speed as measured in the observer 
frame. Here, the expressions of the type [A] stand for the 
difference between the final and the initial value of A. That 
is [A] — Ai — Ao, where the index "0" corresponds to the ini- 
tial state and the index "1" to the final state. For the initial 
state Zo = lp is the pulse width, whereas for the final state 
h = l P + l s is the wavelength of this periodic configuration, 
as measured in the observer frame. 

Denoting the cell mass and magnetic flux as M = pjl 
and m = Bl, the equations of energy and momentum can be 
written as 



M[j] + [pl( K j 2 - 1)] + * 2 [i(l - l/2 7 2 )] = , 



l 



M[yv] + [plKV-y 2 ] + * 2 [y] = . 



(26) 
(27) 



These equations are to be solved for the thermodynamic 
pressure, p%, and the Lorentz factor, 71, of the final state. 
Subtracting them and using the approximation v ~ 1 — 
I/27 2 , one finds 



M 



{2-K)[pl), 



Since po — this yields 



M 



Pi 



(28) 



(29) 



h(2-K 

Next one can simplify Eg 1271 by putting v = 1, which gives 



M[f] + [plu'y 2 ] + * 2 



0. 



(30) 



Elimination of pi from this equation leads to the quadratic 
equation 



fix + (1 — fi)x — (1 + <jq5) = 0, 



(31) 



where x = 71/70, H = T/(2 — T), and 5 = (li — lo)/h = 
5i/(l + 5i). This equation has only one physical solution 



2\x 



(32) 
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ds 2 = (-1 + Po)dt 2 + 2fi dxdt + dx 2 + dy 2 + dz 2 



(37) 



500 



1000 



Figure 9. Energy balance of the adiabatic flow. The dot-dashed 
line shows the fraction of initial magnetic energy turned into heat, 
the dashed line shows the fraction of initial magnetic energy con- 
verted into the bulk motion kinetic energy, and the solid line 
shows the total fraction of utilised magnetic energy. 



where D = (1 — fi) 2 + 4/i(l + ao5). It easy to see that x > 1 
and thus the flow is accelerated. The corresponding increase 
of the bulk kinetic energy is 



, =-(*-!), 
whereas the released thermal energy is 

1), 



[St] m , 

c-m,() <7o 



(33) 



(34) 



both quantities being normalised to the initial magnetic en- 
ergy, £m,0- 

For ero<5 2J> 1, which implies high initial magnetization 
and not very narrow gaps, one has 




» 1. 



(35) 



1/2 

Thus 71 ~ 7oo" , which is significantly lower compared to 
7l = 7of"o, corresponding to the total conversion of magnetic 
energy into the kinetic one. Hence most of the magnetic 
energy must be converted into heat. Indeed, from Eqs. (|33[) 
and p4[) one has 



-[£ k ] > [f k ] 



(36) 



The computational domain covers one wavelength of 
the flow and the initial solution describes a uniform pulse 
located right in the middle of this domain. The gaps are 
not absolutely empty but they are so highly rarefied and 
weakly magnetized that have very little effect on the so- 
lution. The initial Lorentz factor is uniform throughout the 
domain. Both at the left and the right boundaries we impose 
the periodic boundary conditions. 

Figure [8] illustrates the typical evolution of such a flow. 
In this particular model, the initial magnetization of the 
pulse is (Jo = 10, its thermal pressure is negligibly small, 
and it moves with the Lorentz factor 70 = 5 to the right. 
The flow is super fast magnetosonic, with the Mach number 
M ~ 1.6. The pulse (0.25 < x < 0.75) and the gap are equal 
in linear size, with l p — l g = 0.5. The ratio of specific heat 
is again F = 4/3 ( k = 4 and /1 = 2). 

At t = 0.5 the solu t ion is very much as this was antic- 
ipated in iGranot et all l|2010h . Indeed, it is dominated by 
two strong rarefaction waves, both originated at the pulse 
boundaries and moving inside the pulse. However, the one 
produced at the head propagates through the pulse much 
faster than the one produced at the tail. As the result, inside 
most of the pulse the magnetic pressure declines towards the 
head and the magnetic pressure force accelerates the flow. 
The other rarefaction gets stuck at the back where it ejects 
pulse plasma into the tail, which grows in size almost at the 
speed of light. By t = 0.5 the tip of the tail has already 
crossed the gap and collided with head (of the other pulse). 
Two shock waves, one forward and one reversed, are pro- 
duced as result of this collision, they are responsible for the 
spike observed at x ~ 0.77. 

At t = 1 the solution already looks rather different from 
that of an isolated pulse. The reverse shock has moved well 
inside the pulse and it is now located at x ~ 0.63. Behind 
the shock, the flow has almost recovered its initial Lorentz 
factor, magnetization a, and magnetic pressure. 

By t = 10 all the features of the strong rarefaction wave, 
characteristic for the isolated shell solution, have gone. Now 
the solution can be broadly described as a flat-top pulse with 
a tail. 

Due to the periodic boundary conditions, when a wave 
reaches one boundary it re-appears from the other one. The 
forward and reverse waves do this at a very different rate. 
Indeed, if the wave Lorentz factor in the flow frame is 7 W 2> 1 
and the flow Lorentz factor 7 3> 1 then the observed relative 
speed of the forward wave is 



5/3fw = Iftw -/3\~ — 



1 

27 1 



(38) 



whereas for the reverse wave we have 



4.2 Numerical simulations 

In these simulations, we utilise the spatial grid which moves 
relative to the inertial frame of our fiducial observer with 
the initial speed of the pulse. This is similar to the so-called 
"moving window" approach. Namely, we use the time-like 
foliation of space-time defined by the time of observer's in- 
ertial frame, but introduce new spatial coordinate via the 
transformation x — > x — Pot. This leads to the metric form 



5/3rw = |jSrw — P\ 



111 



7 + ll 



(39) 



which is much faster. This difference must be behind the 
observed much more rapid decay of reverse waves compared 
to forward waves. As one can see in Figure [8] at t = 1000 
the amplitude of forward waves is much higher and the flow 
exhibits the characteristic "saw-tooth" profile similar to the 
one which develops in the problem of nonlinear wave steep- 
ening. 
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Figure 10. Differences between the adiabatic flow and the cooling flow with the same initial parameters (<ro = 10, lp = lg 
r = 4/3, and 70 = 5) at the time t = 1000. Solid lines show the cooling flow and dashed lines the adiabatic one. 



0.5, 



Figure [5] illustrates the energy evolution of the flow. It 
shows the variations of the magnetic, bulk motion, and ther- 
mal energies per train wavelength, all being normalised to 
the initial magnetic energy. Since, in this problem S = 0.5, by 
the time the system relaxes to a uniform state its magnetic 
energy decreases by 50 percent. As one can see in Figure 
[5] this has almost been achieved at t = 1000. According to 
the analysis of Sec l4.ll the asymptotic parameters of this 
flow are 7 ~ 270 = 10, a ~ 2.5 and A£k ~ 0.1. These are 
indeed very close to the values observed at t — 1000 (see 
Figures [9] and [8). 

The characteristic time of relaxation towards the uni- 
form state must be given by the dissipation time scale of the 
forward waves associated with the saw-tooth structure of the 
flow. First, the pulse has to pass through the forward fast 
shock where the dissipation occurs. Given the result (|38[) , 
the corresponding time, At c — 2Z7 2 , is actually independent 
on the shock speed, and hence its strength. Second, each 
time when a strong shock crosses the pulse it dissipates only 
a fraction ~ l/2(j of the available pulse energy (see Eq |12l) . 
Thus, the relaxation time can be estimated as 

At T ~ 4laj 2 . (40) 

This turns out to be about the same as the time scale of 
the nonl inear steepening of a fast magnetosonic wave (see 
Eq.ll in lLvuba rskv 2003). For our numerical model, At r ~ 



1000, if we use the parameters of the initial solution for 
the calculations, in excellent agreement with the numerical 
results. 

4.3 Comments on the geometrical effects 

The results of our study show that the interaction between 
individual pulses (magnetic shells) becomes important well 
before the coasting phase of an isolated pulse, significantly 
reducing the efficiency of impulsive acceleration. However, 
the dynamics of adiabatic flows is strongly influenced by 
their geometry. For example, the sideways expansion of con- 
ical jets is an efficient way of converting their thermal en- 
ergy into the kinetic energy of bulk motion. Thus, provided 
the radiative cooling time is long compared to the adiabatic 
one, a conical flow can still eventually become kinetic energy 
dominated. In the context of the pul sar wind accelera tion 
this mechanism has been discussed in lLvubarskvll|2003h . 

In order to find the asymptotic flow parameters in this 
case we simply notice that because 

6n,l = £m,0(l - S) (41) 

the kinetic energy 

£k,l = £k,0 + £m,o5 ~ £ m ,0<5 . (42) 

Thus, the asymptotic magnetization parameter should be 
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Figure 11. Left panel: Energy balance of the adiabatic and the radiatively cooling flows with the same initial parameters. The solid 
lines show the total fraction of utilized magnetic energy (the top line), the fraction of the magnetic energy converted into radiation (the 
middle line), and the fraction of the magnetic energy converted into the bulk motion kinetic energy (the bottom line) for the model 
with fast radiative cooling. The dash-dotted lines show the total fraction of utilized magnetic energy (the top line), the fraction of the 
magnetic energy converted into heat (the middle line), and the fraction of the magnetic energy converted into the bulk motion kinetic 
energy (the bottom line) for the adiabatic model. Right panel: Energy balance of fast cooling flows with different initial Lorentz factors. 
There are three pairs of lines, solid lines for 70 = 3, dashed lines for 70 = 5, and dash-dotted lines for 70 = 10. In each pair, the top 
line shows the total fraction of utilized magnetic energy, and the bottom line shows the fraction converted into the bulk motion kinetic 
energy. The difference between the top and the bottom lines gives the fraction of the magnetic energy converted into the radiation. 



ai ~SpL = Q.-S)/6. (43) 

tk,l 

We comment here that unless the gap is much wider than 
the pulse, and hence S is very close to unity, the asymptotic 
magnetization is still not much lower than unity. For exam- 
ple, when l p = l g , and hence 8 — 1/2, this equation gives 
<7i = 1. Eqs. (|22p and (|23p give us another expression for <ti, 
namely 

<ri =(7o — (1-*). (44) 
7i 

Combining the last two equations we find the asymptotic 
Lorentz factor 

71 = JoCqS . (45) 



5 RADIATIVELY COOLING FLOW 

On the opposite extreme to the adiabatic flow is the case 
where the radiative cooling time scale is much shorter com- 
pared to the time scale of adiabatic cooling. This is the so- 
called "fast radiative cooling" regime, which has been often 
discussed in connection to GRB jets. In this case the dif- 
ference between the plane and spherical geometry is unim- 
portant as all of the released thermal energy is lost to the 
radiation. Assuming that the total fraction of the utilised 
magnetic energy is the same as in the adiabatic regime, one 
would expect very high radiative efficiency, 

7?r = M^. (46 ) 

In fact, this is the same as in the chamber problem with 
empty gaps (see Eq|19[). However, both the radiative cool- 
ing and the radiation reaction force may modify the flow 
dynamics and in this Section we investigate their roles. 



Although in this case we almost recover the asymptotic 
parameters of isolated shells l|Granot et al.|[201ol ) , the accel- 
eration mechanism is different. Indeed, the key role is played 
by the shock heating and thermal acceleration, instead of 
the magnetic pressure acceleration. Moreover, in the con- 
text of GRBs, this regime is not very attractive because it 
implies low radiative efficiency. Indeed, by the time the flow 
becomes kinetic energy dominated, its impulsive origin is 
"washed out", with the remaining shock waves being week 
and allowing dissipation of only a small fraction of the flow 
power. 



5.1 Asymptotic state 

In this case, one can still try to determine the asymptotic 
flow parameters using the conservation laws, just like this 
was done in Sec|3] While the mass and the magnetic flux 
conservation laws remain the same, the laws for energy and 
momentum have to be modified in order to account for the 
radiative losses: 

[( P7 2 + B 2 (1 - 1/2 7 2 ))Z] = Q 7 rad, (47) 
[(P7 2 + B 2 )vl] = Qlrad^rad • (48) 
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Here we use the fact that the total energy-momentum of 
emitted photons can be written as Qu^ad ( see Eq |2UJl , where 
«r ad is the averaged 4-velocity of the flow during its relax- 
ation. Its value depends on possible correlation between the 
rate of radiative losses and the flow speed, e.g. the radiation 
may come mostly from the fastest parts of the flow. Since 
this is essentially an unknown parameter, we end up with 
an under-determined system, which has only four equations 
and five unknowns, namely 71, Bi, pi, 7 ra( j, and Q. 

Eliminating Q, B\, pi, and using the usual high speed 
approximation v ~ 1 — I/27 2 , we find the equation deter- 
mining 71 as a function of 7 ra d, 

teX-f )-(-$)—• « 

If we put 7rad = 71 then this equation yields 

71 = 70(1 + <7o<5) 1//2 ~ 70V aoS , (50) 

where the last step assumes aoS 3> 1. This is indeed very 
similar to what we have found in the adiabatic case in the 
slab geometry (see Eq |35|) . 

On the other hand, when 7 ra d = 70, we have 

71 = y(<ro<5 + {°l8 2 + 4) 1/2 ) ~ j a S . (51) 

This implies very efficient flow acceleration with almost total 
conversion of magnetic energy into the bulk motion kinetic 
energy for S ~ 1, contrary to what we have anticipated. The 
reason for this is the strong radiation reaction force, which 
accelerates the flow. Indeed, 7 ra d = 70 implies that most of 
the time 7 ra( j is lower compared to the center-of-momentum 
Lorentz factor. Because of this, in the centre-of-momentum 
frame, the photons are emitted mostly in the direction oppo- 
site to the flow direction. This seems hardly possible as most 
of the radiation must come from the immediate downstream 
of forward shocks and if the flow preserves the structure of 
adiabatic solution then these are the fastest sections of the 
flow (see Figure [8]). In fact, this argument suggests that the 



asymptotic Lorentz factor can be quite close to that given 
by Eq |5UI In order to verify this deduction we have carried 
the numerical simulations described below. 

5.2 Numerical simulations 

In these simulations we used the same cooling function 
as in the chamber problem (see Eq|21[l, In order for the 
numerical shock structure not to have much influence on 
the overall effect of radiative cooling, we had to insure 
that the cooling length scale was significantly higher com- 
pared to the shock thickness. If the proper cooling time 
scale is Ar CO oi then in the observer frame this scale is 
At coo i = 7Atcooi. In this time, a weak forward magnetosonic 
shock wave moves relative to the flow by the observed dis- 
tance AZ ~ (l/7 2 )At coo i ~ At coo i/7. Thus in order for the 
length of the cooling region behind the shock to be a frac- 
tion a = AZ/Z of the flow wavelength, the proper cooling 
time should be Ar CO oi = «7Z. In these simulations we used 
a = 0.2. 

First we studied the model with the same initial param- 
eters as the adiabatic one described in Section [4] ao — 10, 
l p = l g — 0.5, r = 4/3, and 70 = 5. Overall, the cooling 
flow shows a similar deviation from the d ynamics of isolated 
pulse, considered in iGranot et ail (|2010n . as the adiabatic 
model. Figure [TU] compares the adiabatic and cooling flows 
at later times. As one can see, both flows develop the char- 
acteristic "saw-tooth" profiles and apart from the tempera- 
ture, their parameters are rather similar. The cooling flow 
is only slightly faster compared to the adiabatic one. 

The left panel of Figure [11] compares the energy balance 
of both flows. It confirms that the cooling flow is slightly 
more efficient in converting the magnetic energy into the 
bulk kinetic energy. It also shows that by t — 1000 most of 
the free magnetic energy has been utilized. Comparing the 
numerical results with the predicted asymptotic parameters 
of the cooling flow (see Sec l5.1fl , we find that Eq |50l does 
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Figure 13. Evolution of the Lorentz factor in the model with 
filled gaps. The curves show the solutions at t = 1 (dotted line), 
t = 100 (dash-dotted line), t = 1000 (dashed line), and t = 3000 
(solid line). 



better, with 71 ~ 11, whereas Ea l51l significantly overesti- 
mates the asymptotic Lorentz factor, giving 71 ~ 50. Figure 
[10] shows that most of the radiative energy losses are asso- 
ciated with the fastest parts of the flow, just like we have 
anticipated, and this explains why Eq |50l provides a more 
accurate estimate. 

The right panel of Figure [11] compares the energy bud- 
gets of three cooling models which differ only by the Lorentz 
factor of the initial solution, 70 = 3, 5, and 10. One can see 
that the efficiency of magnetic acceleration grows with 70, 
but only slightly. In all these models, more than half of the 
free magnetic energy is converted into radiation. Supported 
by the strong arguments presented in Sec l5.ll we conclude 
that when oq8 3> 1 most of the released magnetic energy 
is converted into radiation, and the radiative efficiency can 
indeed be estimated using Eg 1461 



Here we present the results for the model with param- 
eters 70 = 5, <T P = 10, <7g = 0.1, Si = 5 m = 1, and 8 C — 0.01. 
For this parameters Eg 1521 gives [£ m ] = — 0.4£ m .o. The left 
panel of Fig |12l compares the energy balance of this model 
with the empty gap model which has the same pulse param- 
eters and hence [£ m ] = — 0.5£ m ,o. As one can see, in both 
these models the released magnetic energy does evolve to- 
wards the predicted asymptotic values and at approximately 
the same rate. The partition of this energy between the bulk 
kinetic energy and the radiation is also similar. This allows 
us to conclude that in the case of energetically subdominant 
gaps, the radiative efficiency is similar to what we have for 
empty gaps, and is well described by 



V 



(53) 



(l+(5l)(l+(5e) ■ 

Because the total accelerated mass is twice as higher in 
the model with filled gaps, one would expect the asymptotic 
Lorentz factor in this model to be lower compared the cor- 
responding model with empty gaps. The data presented in 
Fig |13l and Fig llOl show that this is indeed the case. How- 
ever, the difference is rather small, and this suggests that 
the asymptotic Lorentz factor is still well described by Eo l50l 
when the gap mass does not exceed that of the pulse. 

The energy curves shown in the left panel of Fig |12l are 
more ragged for the model with filled gaps. Moreover, their 
structure suggests a quasi-periodic process. In fact, the curve 
of radiated energy is similar in shape to those found in the 
chamber problem, indicating a strong variation of the dis- 
sipation rate. This is confirmed by the left panel of Fig |12l 
which shows the energy loss rate. Just like in the cham- 
ber problem, one would expect its spikes to be associated 
with crossings of weakly magnetised gaps by shocks. This 
is indeed the case, as illustrated in Fig |14l This plot also 
shows that, like in other models, the flow develops the char- 
acteristic saw-tooth profile. The time separation between 
these spikes can be explained using Eo l38l According to 
it, the pulse crossing time by a relativistic forward wave 
is At c ~ 2j 2 l ~ 50, where we used 7 — 70 = 5. This is in- 
deed very close to the spike separation in Fig |12l The rapid 
variability inside the first spike is connected to the reverse 
shock, which crosses the shell much faster (see Eq |39[) . 



6 RADIATIVELY COOLING FLOW WITH 
FILLED GAPS 

Finally, we investigate a flow with fast radiative cooling and 
gaps filled with low magnetised plasma from the start. Pro- 
vided the initial state has the same Lorentz factor both for 
the gap and the pulse, the fraction of released magnetic en- 
ergy is the same as in the chamber problem (see Eq |18|l 

M = 1 d + ^ 1/2 ) 2 l52) 

£m,0 (l+&)(l + ftO ■ V ' 

This energy is converted partly into the radiation and partly 
into the bulk kinetic energy. Although this partition is dif- 
ficult to estimate analytically, there are no obvious factors 
that could significantly shift the balance in favour of the ki- 
netic energy, suggesting that most of the magnetic energy 
should still be converted into the radiation. This is confirmed 
by numerical simulations 



7 DISCUSSION 

7.1 Dynamics of impulsive 

magnetically-dominated outflows 

In recent years a great deal of progress has been achieved in 
the dynamics of magnetic steady-state flows but the observa- 
tions seem to suggest that the central engines of cosmic jets 
are highly variable. Is this variability only a minor compli- 
cation or it can actually lead to a qualitatively different jet 
dynamics? The rece nt pioneering studies of this issue have 
suggested the latter (IContopouloslll995l;lG"ranot et alj|20ld : 
ILvutikovl |2010| ; iLvutikov fc Listed |2010| ; iLvutikovl 120111 ). 
However, they were mainly concerned with the dynamics 
of individual shells, presumably ejected by the central en- 
gine, and their interaction with the external medium. These 
issues were studying using rigorous mathematical modelling. 
On the contrary, the problem of interaction between mag- 
netically accelerated shells has remained until now a subject 



© 0000 RAS, MNRAS 000, 000-000 



Shock dissipation 17 




Figure 14. Solution for the model with filled gap at t = 80, 120, and 180. The solid line shows the rest mass density and the dashed 
line shows the total pressure, dominated by the magnetic pressure. As one can see, at t = 80 and 180 the fast forward shock crosses the 
flow section which has significantly higher rest mass density, the gap section. At the same times, the radiation loss rate, which is shown 
in right panel of Figure [T2l exhibits strong peaks. 



of only rather speculative semi-quantitative analysis. One of 
the main objectives of our study was to explore this impor- 
tant issue a bit further. 

One of the most important properties of the isolated 
shell solution is the concentration of energy and momentum 
in the compact head o f the shell formed b y the reverse rar- 
efaction wave early on (|Granot et alj|20ich . This rarefaction 
wave can be seen in the plots of the left column of Figure [HI 
where it occupies the domain 0.45 < x < 0.75. The head for- 
mation is completed when this rarefaction reaches the back 
of the shell and gets reflected as a forward rarefaction. The 
part of the shell located between the leading interface with 
vacuum and the leading front of the reflected rarefaction 
wave is what constitutes the shell head. In the source frame, 
it shows very little spreading in the direction of motion and 
its mean Lorentz factor increases as r 1 ^ 3 . 

iGranot et all (|20ld ) and iGranotl l|2011br i argued that 
this somewhat unusual, but well understood in the relativis- 
tic framework, behaviour of the shell head allows us to ignore 
the interaction between individual shells until their Poynting 
flux is almost fully converted into the bulk motion kinetic 
energy. Only after this, during the coasting phase, the shell 
spreading becomes important and leads to strong collisions 
between shells, thus allowing internal shocks with high dis- 
sipation efficiency. This assertion had apparently given 
a great importance to the assumption, made in that 
study, that the gaps between shells were empty, or 
at least so highly rarefied that their plasma could 
not influence the shell dynamics. In fact, the au- 
thors overlooked the ejection of shell plasma into its 
tail, which would rapidly fills the gap, and did not 
explore the possible implications of the head-tail col- 
lision for the shell dynamics. 

An isolated shell enters the coasting phase at t c ~ 
(£p/c) 7^ where 7 e is th e asymptotic Lorentz factor of the 
shell l|Granot et al.ll2010h . This is the time required for a fast 
magnetosonic wave moving relative to the shell head with 
the observed speed 5v^ = c/j% (see Eq |38|) to traverse the 
head. The relative observed speed at which the shell plasma 
is ejected into the tail can be estimated via Eq |39l which 



gives us 5vt ~ 2a(j c /y) 2 5vh ~Ss> Svh- Thus, one can avoid the 
head-tail collision only if the gaps between shells are much 
wider than the shells. 

Our simulations have shown that, at least in the 
case there shells and gaps have comparable widths, 
the tail-head collision strongly modifies the flow dy- 
namics, making the results for an isolated shell ir- 
relevant. The gaps soon become filled up with sig- 
nificant amounts of plasma, and thereafter we have 
what is best described as a continuous inhomoge- 
neous flow superimposed with a wave train. A sim- 
ilar outcome is observed when the gaps are filled 
with plasma already from the start. The subsequent 
acceleration and radiation of this flow is determined 
by the dissipation rate of shocks, which form an in- 
tegral part of the train, the plasma cooling rates, 
and the flow geometry. In any case, the fraction of 
magnetic energy that can be converted into either 
the energy of radiation or the kinetic energy of the 
flow is set by the energy of the train. Under the con- 
dition of magnetic flux freezing, it is simply given by 
the degree of expansion in the longitudinal direction 
which the magnetised shells can achieve during the 
transition to force equilibrium and thus by the am- 
plitude of initial fluctuations of magnetic pressure. 

When the radiative cooling can be ignored, the asymp- 
totic flow parameters follow from the basic conservation laws 
of mass, magnetic flux, energy, and momentum. This way 
we find that in the case of slab geometry, the asymptotic 
Lorentz factor is 7 ~ 7o(o"o<5) , rather than 7 ~ 70 uo 
expected in the case of the total conversion of the mag- 
netic energy into the bulk motion kinetic energy (Here 
5 — l g / (Z g + l p ) is the relative thicknes s of the gap .). Incon- 
trast to the case of an isolated pulse IGranot et "all l|2010h . 
the magnetic acceleration is inefficient and most of the free 
magnetic energy is converted into the thermal energy via the 
shock dissipation. However in the case of spherical geometry, 
the sideways expansion of the flow is accompanied by effi- 
cient conversion of this thermal energy into the bulk kinetic 
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energy, via the thermal acceleration mechanism. This leads 
to the asymptotic Lorentz factor 7 ~ 7o<to<5. 

When the radiative cooling is much faster than the adi- 
abatic one, the difference between flows with slab and spher- 
ical geometry is unlikely to be significant, although the pres- 
ence of poorly radiating components such as protons com- 
plicates this issue and requires further investigation. Even 
without such components, the fast cooling regime is more 
involved as the conservation laws alone do not determine 
the asymptotic flow parameters. The outcome now depends 
on details of the radiative cooling process. The radiation re- 
action force may both help the flow acceleration and make it 
more difficult, depending on whether the photons are emit- 
ted predominantly in the forward direction in the centre-of- 
momentum frame or in the backward direction. From the 
numerical simulations we find that the emission comes from 
the fastest part of the flow and thus the radiation reaction 
force is a decelerating one. Most of the utilised magnetic 
energy is converted into the radiation and the asymptotic 
Lorentz factor of the flow remains significantly below the 
value of 7oo"o5, characteristic of efficient magnetic accelera- 
tion. 



7.2 Astrophysical implications 

In terms of astrophysical implications, our results are most 
relevant to the issue of internal shock dissipation and emis- 
sion from relativistic jets of GRBs and AGN. A detailed test 
of the shock dissipation model in magnetically-dominated 
flows against observations is beyond the scope of this paper. 
Here we only outline few issues to be investigated in future 
studies. 

Just like in the hydrodynamic model, the energy reser- 
voir for the internal shock dissipation in the magnetic model 
is associated with the variable component of the jet. For 
Poynting-dominated jets this is the magnetic energy of the 
fast magnetosonic waves driven into the jet by the variable 
central engine. Their contribution to the overall jet energy 
budget depends on the exact details of the central engine op- 
eration, and at least in principal it may well be dominant. 

Provided the radiative cooling time is sufficiently short, 
most of the dissipated energy is converted into radiation. In 
our simplified test problems, the radiative efficiency is given 
by the Equation 1181 which shows that it can be very sub- 
stantial even for a rather moderate central engine variabil- 
ity. Thus as far as the observed radiative efficiency of GRB 
jets is concerned, the internal shock model of prompt GRB 
emission and the magnetically-dominated model of GRB jets 
are not mutually excluding. However, the dissipation length 
scale increases with the jet magnetization and may exceed 
the jet length. 

In the hydrodynamical version of the internal shock 
model, the prompt emission originates from the region where 
individual ballistic shells collide wit h each othe r because of 
the differences in the ejection speed (|Piranll2004h . For strong 
variation of the Lorentz factor, this occurs at 

R B ~ yfcSts ~ 3 x 10 16 cm (^j , (54) 

where St s is the characteristic time interval between the shell 



ejection^], and jj is the Lorentz factor of the slower shell. 
Each collision gives rise to a pulse of the prompt emission 
light-curve. Its duration, St p , is determined by the curvature 
of the shock front and the Doppler beaming (the so-called 
"angul ar spreading effect" ) , and it turns out to be the same 
a.s 5t e (|Piranll2004l ). 

In the magnetic model, the prompt emitting region can 
extend well beyond Rs. According to Eq(40]a forward shock 
dissipates its energy when it covers the distance R ~ aRa. 
For t p ~ 1 s, R s is already dangerously close to the radius of 
the external shock which seems to rule out a 3> 1. However, 
the problem arises only if we associate each strong individ- 
ual pulse of a GRB light curve with an individual shock 
wave. If instead each such pulse is associated with a whole 
packet of shocks, so that 8t s < St p /a, this is no longer an 
issue. Obviously, this explanation implies a secondary cen- 
tral engine variability process, operating on the time St p , 
which modulates the output of the primary process, oper- 
ating on the timescale St s . In fact, this way one can explain 
why about 2 0% of GRBs have rather featureless light curves 
(|Piranll2004h . This may just be the case of weak modulation. 

In this model, there are ~ a shocks inside the dissi- 
pation zone and one may wonder if their emission signals 
overlap. In the source frame the relative speed of a light sig- 
nal and a forward fast magnetosonic wave propagating in 
the radial direction is 

*M * ^2 • (55) 

The travel time across the dissipation zone is aR s /c. During 
this time the distance between these two waves is changed 
by only c5t s /8. Thus, the overlapping could only be a result 
of the angular spreading effect, which spreads 5-shape sig- 
nals over the time 5t ang = R/2cy^ , where R is the emission 
radius. This gives us <5i ang = 5t s /2 at the beginning of the 
dissipation zone (R ~ R s ) and <5i an g = aSt a /2 3> 5t s at the 
end of the dissipation zone (R ~ <jR s ). One may draw two 
conclusions from these numbers. First, a micro-pulse from 
an individual shock must have a long smooth tail of length 
~ a5t s /2. Spike-like features on the scale of ~ St a can be 
found only at the head of this micro-pulse (They may cor- 
respond to shock crossings of contact discontinuities in the 
jet; see Fig(§] and FigfTT}. Second, the tails of individual 
micro-pulses overlap but their leading spikes do not and can 
give rise to a pulse substructure. The fact that, shocks are 
strongest at the beginning of the dissipation zone increases 
chances of them being detected. This could be the origin of 
the observed variability o f GRBs on millisecond time scale 
(| Walker fc Schaefer|[200ol ). 

How short can 5t a be? In the Blandford-Znajek model 
of GRB central engine, the only "easy" way of changing 
its jet power, Lj , is via changing the flux of open magnetic 
field lines, "J, which threa ds the black hole, as these p aram- 
eters relate via Lj oc * 2 |Blandford fc Znaiek|[l977f) . This 
magnetic flux may change significantly if the disk drags in 
an alternating magnetic field, which may even bring about 
a change of polarity of the black hole magnetospherfl The 

1 Such a characteristic time does not have to exist. Instead, the 
central engine may exhibit a wide distribution for the shell ejec- 
tion time. 

2 The effective loss of magnetosphere and its shielding action 
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relevant time scale for this process is probably the disk inner 
edge "viscous" time scale. For an a-disk, this is 



St v « 10 



a5 2 
10" : 



M 



(56) 



where M is the black hole m ass and 5 = Hd/Rd is the r atio of 
the disk height to its radius (jShakura fc Sunvaevll 973'). This 
is almost a hundred times shorter comp ared to the observe d 
typical separation between GRB pulses dNorris et al.lll996h . 
horizon. The shortest scale for restructuring of the black 
hole magnetosphere is given by the light crossing time of 
the ergosphere 



ft. = 2 ™ „ W 



M 

Mo 



/ts 



(57) 



A rapid and frequent restructuring of magnetosphere may 
also be typical for a newly-born millisecond magnetar, due to 
its ultra-strong magnetic field, lack of solid crust, and active 
magnetic dynamo. These estimates assume fast magnetic 
reconnection in the magnetospheres, which has been ques- 
tioned recently on the basis of collisi onal nature of reconnec- 
tion o f super-strong magnetic field jMcKinnev fe Uzdenskvl 
l20ld ; iLvutikov fc McKinnevI 1201 lh lUzdenskvl 1201 ll ). al 
though even the collisional magnetic reconnection can be 
fast due to the second ary tearing insta bility (see the discus- 
sion and references in IUzdenskvll20Tlh . 

As to the modulation process, it should also involve 
variation of <]/ but on a larger time scale. For example, more 
massive accretion disks could support stronger magnetic 
field and result in larger magnetic flux trapped by the black 
hole. Variations of the disk mass could result from unsteady 
stellar collapse. The typical duration of strong GRB pulses 
is around the free-fall time from the radius of ~ 10 9 cm for 
the GRB progenitor in the collapsar model. Given the total 
radius of the progenitor, ~ Rq, it does not seem implausible 
for the GRB jet to perturb the stellar mass distribution on 
this scale. 

Swift observations of early X-ray afterglows have dis- 
covered the presence of plateaus as well as strong flares in 
their light-cu rves, both une xpected in the original external 
shock model (|ZhangH2007br ). One possible and perhaps the 
most likely explanation of these features is that they rep- 
resent t he "late prompt emis s ion" of a long-lasting central 
engine (|Ghisellini et al.ll2007l . I2009T ). This interpretation is 
supported by a number of sim ilarities between the X-ray 
flares and the GRB sub-pulses (jMargutti et al.ll2010T ). Such 
a long-lasting activity, at least up to 10 4 seconds, is incon- 
sistent with the high mass accretion rate, above lO -2 M0/s, 
required by the neutrino annihilation mechanism of jet pro- 
duction and implies the magnetic mechanism. 

In addition to the similarities between the X-ray flares 
and the GRB sub-p ulses, there are also few differences 
jMargutti et al.ll2010h . For example, the flare duration in- 
creases linearly with the flare time since GRB trigger. In 
our model, the duration of prompt pulses is determined by 
the duration of strong accretion episodes in the history of 
the central engine. If the X-ray flares are associated with the 



during a change of polarity may also let weakly magnetized sur- 
rounding plasma to enter the jet channel and to become entrained 
by the jet. 



fallback accretion, then the flare time is likely to be deter- 
mined by the location of fallback turning point whereas its 
duration by the spatial dispersion of falling back material. 
They may well correlate with each other. 

The evolution of individual pulses of GRBs is of- 
ten described a s "Fast Rise Exponential Decay" (FRED, 
e.g. |Piranll2004h . The latest results show that on average 
both the prompt pulses and the X-ray flares have approx- 
imately twi ce as shorter the rise time compared to the 
decay time (|Norris et al.ll2005l ; Ichincarini etahl |2010| 1 ■ Al- 
though shorter rise times are expected in the standard in- 
ternal shock mode l , it d oes not really say by how much 
l|Yi-Ping Qin et all 120041 ') . Perhaps, this tells us that the 
pulse shape is not determined by the shell collision after 
all. Other processes may have similar time properties. In 
fact, FRED mass accretion rate could be a general conse- 
quence of the diffu sive transport in episodic accretion disks 
i|Wood et al.ll200 il l. 

We have shown that, provided the GRB jets are cooled 
mainly radiatively, they are likely to remain Poynting- 
dominated even beyond the prompt emission zone. The ex- 
act nature of piston driving the external shock into the sur- 
roundings of GRB has little effect on the shock dynamics 
and emission. Provided the total energetics is the same, both 
the kinetic energy dominated and the magnetic energy domi- 
nated ejecta pro duce the same afterglow emission associated 
with t his shock (ILvutikov fc Blandfordf2003l ; lLvutikovl20"05l 



l201ll ; iMimica et al. I 120091 '). However, the reverse shock, 



which is driven into the ejecta, will b e much less dissipa- 
tive if the ejecta is highly magnetized II Zhang fc Kobavashil 
l2005l;lMimica et al. |[2009h . This could be the reason behind 
the paucity of the lower energy, most likely optical, flashes 
expected from the reverse shock in the non-magnetic models 
of afterglows l|Gomboc et al.ll2009h . 

As to the issue of the jet acceleration, our results in- 
dicate that the impulsive magne tic acceleration mechanism, 
proposed bv lGranot et al.l l|2010l ). is unlikely to operate in its 
original form. Even if we ignore the entrainment of surround- 
ing plasma, the strong shock interaction between heads and 
tails of magnetic shells, revealed in the present study, is al- 
ready sufficient to prevent the conversion of almost all avail- 
able magnetic energy into the kinetic energy shell heads. In- 
stead, the problem reduces to its more standard form where 
non-linear waves travel along the jet, interacting with the 
mean flow. The efficiency of the mean flow acceleration by 
these waves depends on the efficiency of its radiative cooling, 
and hence on the photon opacity of the flow. If this cooling 
is weak then the jet acceleration can be very efficient, with 
a large fraction of the wave Poynting flux converted into 
the kinetic energy of the flow. Otherwise, it is converted 
into radiation. Based only on baryonic electrons the optical 
thickness of a GRB jet to Thomson scattering is 



L 6 



ri3<T7j,2 



(58) 



where L is the isotropic jet power. Thus, it seems unlikely 
for the dissipation zone to be optically thick, unless pairs 
dominate the opacity. This agrees with the high observed 
radiative efficiency of GRB jets. 

Impulsive jet production provides favourable conditions 
for entrainment of weakly magnetised plasma that may ex- 
ist in the vicinity of the central engine. We focused on the 
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implications of this process for the dissipation efficiency of 
shock waves in magnetically dominated flows. But in addi- 
tion to this, the presence of weakly magnetised plasma may 
also help to overcome another difficulty of the shock model, 
related to the shock acceleration of non-thermal particles. 
However, for this to work the magnetization may have t o be 
very low, down to a < fO -3 (jSironi fc Spitkovsky||201oh . 



8 CONCLUSIONS 

In this paper we analysed the potential role of shock dissi- 
pation on the dynamics and emission of impulsive Poynting- 
dominated relativistic jets. The main insights came from an- 
alytical and numerical solutions of Relativistic MHD equa- 
tions in slab geometry. For computational reasons the nu- 
merical simulations were limited to flows with Lorentz fac- 
tors which were much lower compared to those deduced from 
the observations of GRBs. These and other limitations of 
this study warn against making firm conclusions and a cer- 
tain degree of uncertainty definitely remains in many re- 
spects, and particularly when it comes to astrophysical ap- 
plications. Keeping this in mind, our main conclusions are 

• The dissipation efficiency of strong shocks in highly 
magnetized plasma is low ~ 0.5(1 + er) -1 , mainly because 
only the kinetic energy dissipates and it represents only a 
small fraction of the total energy flowing through the shock. 
For moderate magnetization, a ~ I , the shock dissipation 
efficiency is reasonably high, ~ 30% of the total energy flux. 

• The dynamics of a magnetic shell in a train of shells, 
all ejected in the same direction and initially separated by 
empty space, can be rather diffe rent from that of a n iso- 
lated shell in vacuum, studied in iGranot et all (|2010h . Un- 
less the separation between these shells is very large, they 
strongly interact with each other. The result of this inter- 
action is best described as an underlying continuous flow 
superimposed with strong magnetosonic waves travelling in 
the same direction. The wave Poynting flux provides energy 
for heating and acceleration of the flow. A similar outcome 
is observed in models where gaps are filled with energetically 
sub-dominant plasma from the very beginning. 

• For an infinite flow, the reservoir of magnetic energy 
is given simply by the decrease of the magnetic energy, as 
dictated by the magnetic flux freezing, during the transition 
to a wave-free final state. The shock dissipation is an essen- 
tial part of this transition. The radiative cooling allows to 
maximise the released magnetic energy. For radially diverg- 
ing flows, the adiabatic cooling is expected to do the same. 
In both cases, the fraction of released magnetic energy does 
not depend on the initial flow magnetization but only on 
the initial magnetic inhomogeneity of the flow, or in other 
words on the fraction of energy carried by the waves. How- 
ever, the flow magnetization determines the tempo and the 
characteristic length scale of the energy release, which grows 
oc a for ct 3> I. 

• With application to GRB jets, this shows that the inter- 
nal shocks may still be responsible for the prompt emission 
even in the case of Poynting-dominated jets. At least, the 
high observed fraction of the prompt emission in the total 
energy budget of GRB events is not inconsistent with this 
model. 



• However, the increased size of the prompt emission (in- 
ternal shock dissipation) zone, compared to that of the non- 
magnetic model, calls for a somewhat different interpreta- 
tion of the prompt emission variability. In particular, the 
observed individual pulses with duration around one second 
are unlikely to be associated with the emission of individ- 
ual shocks, as this would push the dissipation zone beyond 
the external shock radius. Instead, each such relatively long 
pulse could represent the combined emission from a whole 
pack of shocks associated with mini-shells ejected by the cen- 
tral engine on a much shorter timescale. The long timescale 
variability would then be related to some modulation pro- 
cess, determining the energy of emitted mini-shells. In fact, 
both the magnetar and black hole magnetospheres do al- 
low variability on millisecond, and possibly even shorter, 
timescales. As to the nature of the modulation, one possi- 
ble mechanism is the unsteady mass supply to the accretion 
disk in the collapsar scenario. 

• The emission from individual shocks may produce fine 
substructure, on the timescale of shell ejection. This could 
be the origin of the observed millisecond spikes in the light- 
curves of the prompt emission. 

• After leaving the zone of internal dissipation, the ini- 
tially Poynting-dominated jet is most likely to remain 
Poynting-dominated. This may not effect the dynamics of 
the forward shock and its afterglow emission, but will result 
in reduced emission from the reverse shock compared to the 
non-magnetic model. Yet, the presence of weakly magne- 
tized domains in the jet will increase the radiative efficiency 
of this shock compared to the value expected for a uniform 
jet with the same mean magnetization. 
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APPENDIX A: RELATIVISTIC 
PERPENDICULAR MHD SHOCKS 

Here we analyse MHD shock waves in the special case where 
the magnetic field is parallel to the shock front and the flow 
velocity is perpendicular to it. In the shock frame, the fluxes 
of energy, momentum, rest mass, and magnetic field are con- 
tinuous across the shock 

(w + B 2 )-y 2 v = const, (AI ) 

B 2 

(w + B 2 )-y 2 v 2 +p+ — = const, (A2) 

p"/v — const, (A3) 

B'yv = const, (A4) 

where p is the rest mass density, p is the gas pressure, w = 
p+ up is the relativistic enthalpy, ft = r/(r — I), where V is 
the adiabatic index, B is the magnetic field as measured in 
the fluid frame, and 7 is the Lorentz factor. We select the 
frame where the velocity vector is normal to the shock plane 
and the magnetic field is parallel to it. We use subscripts I 



© 0000 RAS, MNRAS 000, 000-000 



Shock dissipation 21 



and 2 to denote the upstream an the downstream states 
respectively. 

Equations (|XTj) , (|A3[l and (E4)l yield 

pi _ B'2 _ U'2 _ S 

Pi Bi G\ x 



0% = ~ [ 6(1 + KCtl + <7l) — <7l I — I — 1 



(A5) 
(A6) 



(A7) 

(A8) 

(A9) 
(A10) 



6_ 

« V ' \x 

where 
a = B 2 /p 

is the magnetization parameter, 
a 2 =p/p 

is the temperature parameter, 
X = Wt/vi, 

= + „?(!- 

72 

and u = W7. Using the definitions of a and a, Equation (|A2jl 
can be written as 



(1 + ko ± + (Ji ) h pi (a! + -<7i) = 

pi 2 

7 2 1 
(l + Kal + 0-2)— +p 2 (ai + -^2), (All) 
P'2 2 

where j = py/3. Via substituting the expressions for pi, 02, 
and a,2 from Eqs. (|A5IA6|l this equation becomes an equation 
for Xi which defines it as a function of the upstream state 
parameters. In general, this is a rather combersome algebraic 
equation which has to be solved numerically. However, in 
many astrophysical applications one may assume that the 
upstream state is cold (a\ — \ 0), which allows significant 
simplifications. Then Eq lAl II reduces to 



a-ix +0-2X + a ix + a o +g(x) = ■ 

where 



as = (1 + ai)ui 



(A12) 



9 fM 
«3X + (as + a 2 )x + («3 + a 2 + ai) = 



(A14) 



where 



f(x) = 



g(i)-g(x) = 1 - xS(x)) 
x-i ' (i-x) ' 



For this trivial solution both the numerator and the denom- 
inator in the above expression for f(x) vanish. An addi- 
tional continuous solution, x = 1) mav exist when the shock 
speed relative to the upstream state equals to the fast mag- 
netosonic speed of this state. In order to verify this, one can 
replace /(x) m Eq |A14l with its limiting value 

lim/(x) = l-u\ 

and then substitute x — 1 into the left-hand side of this 
equation. The result is 



Uj = 01, 



(A15) 



which is the fast magnetosonic condition. Indeed, in the limit 
of cold flow, the fast magnetosonic speed is isotropic and 
equals to the Alfven speed 



2 B 

C{ = 



B 2 + p 1+a 



(A16) 



The corresponding Lorentz factor 71 = \/l + a. Thus the 
condition (|A15|) becomes the condition 



Mi = 1 



(A17) 



on the shock fast magnetosonic Mach number, which is de- 
fined as 



Mi = 



Wf,l 



(A18) 



where m = Cfjf. With Mi and o\ as parameters, Eq |A14l 
reads 



ax +b X + c = 2f( X ), 
where 

a = 2(l + ffi)ffi(K - l)Mj 2 



(A19) 



0,2 — — (T1M1 



2k 



ft j \ - 1 ai 
,-(l + ffi)tti - — 



Oi = 



(l + 0-i)(l+M 2 ) 



fflQ 



K-2 

2k 



(1+«i), 



and g(x) — — <^(x)x/ K - O ne can see that the solutions of 
this equation are parametrised only by the upstream mag- 
netization and the shock speed. Instead of the shock speed 
one can introduce the more traditional shock Mach number, 
that will be done later. 

For any values of upstream parameters Eq |A12l must 
allow the trivial continuous solution, x = lj which describes 
a flow without a shock. Thus, we have 

a 3 (x 3 - 1) + o 2 (x 2 - 1) + 01 (x - 1) = 3(1) ~ 9(X), (A13) 
which can also be written as 



b — — <7i(2 + kcti)Mi — wi, 

C = -<7 2 (K - 2)M 2 - KG 1 - 2((Jl + 1). 

For Mi ^> max(l, 1/cri, l/y / oT) one can only retain the 
terms that are proportional to M 2 in the coefficients of x'- 
This yields 

2(l + a l )(K-l)x 2 -(2 + Ka 1 ) X + (2-K)a 1 = -A/( X ).(A20) 

The factor Mj -2 on the right-hand-side term is small, sug- 
gesting that this term can be also ignored. This gives us 
simple quadratic equation which does not involve Mi, 

2(1 + <ti)(« - l) x 2 - (2 + «o-i)x + (2 - k)<ti = 0. (A21) 

Since, the physical meaning of x dictates < x ^ li the 
only suitable solution of this equation is 



X = 



2 + K(Ji + \[~D 
4(k- l)(l + o-i)' 



(A22) 
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where 



D = (2 + ko-i) 2 - 8(1 + (Ti)(Ti(k - 1)(2 - 



T his is the same solution as in the analysis of oblique shocks 
in lKomissarov fe Lvutikov I i|201ll ) and, f or k = 4 (T = 4/3), 
in the perpendicular shock solution of iKennel fc Cororiiti 
(|l984h . 

Further simplification is possible when <ri ^> 1. In this 

case 



1 + 



2k 



3k 



r°"i 



which for k — 4 gives us 
Using this result one can find 

1/2 

72 ^ , 
02 ^ Ml/3l, 



?^2 

B 2 ^ Mi Si 



p m ,2 — MiPnj,!, 

(T2 ~ MlCTl. 



(A23) 

(A24) 

(A25) 
(A26) 

(A27) 

(A28) 
(A29) 

(A30) 

(A31) 



These allow to find how energy is distributed between its 
different forms in the downstream flow. Denoting as = 
fr/ 2 v, Fjn — B 2r y 2 v, and F t = tqrfv the fluxes of kinetic, 
magnetic, and thermal energy respectively, we have 

1 



F\ 



Mi 



1 



Ft,2 

and 

Fm 2 — Jm,l 



(A32) 
(A33) 

(A34) 



Thus, approximately one half of the upstream kinetic energy 
is dissipated into heat, whereas the other half is converted 
into the magnetic energy. The shock dissipation efficiency 
can be denned as 

(A35) 

where Ftat = Fm + Fk + Ft is the total energy flux. From 
the above results it follows that 
1 



Vs 



2(1 + ai ) 



(A36) 



Thus, the shock dissipation efficiency is greatly reduced in 
highly magnetised plasma. 

When (7i > 1 we have f(x) — — 2<7iMl. Substituting 
this result into Eq |A20l one can verify that the right hand 
side term of this equation is indeed much smaller compared 
to any term on the left hand side. 

We should point out th at the condition m 3> 1, stated 
in lKennel fc Coronitil ll 19841 ') for their approximate shock so- 
lution, is not quite correct. The proper condition which leads 
to the solution (IA22|) is Mi 3> 1 and it requires tti *Jai. 



The shock solution does not even exist when Mi < 1 whereas 
the condition ui ^> 1 ma y still be satisfied. 

IZhang fc Kobavashil (|2005h use the Lorentz factor of 
the relative motion between the states on both sides of the 
shock, 

712 = 7172(1 — V1V2), 

as a shock strength parameter. Using the above results, it is 
easy to show that for Mi 3> 1 and <ti ^> 1 

Mi 



712 



and hence this parameter is equivalent to the shock Mach 
number. However, for lower Mi and 01 the connection be- 
tween these parameters is more involved. The shock Mach 
number is a traditional parameter which is rightly recognised 
as most useful in characterising the shock strength. It de- 
scribes the state upstream of the shock, leaving all the down- 
stream state parameters to be found from the shock equa- 
tions, in contrast to 712, which involves both these states. 
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